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A theoretical study of free space coupling to high-Q whispering gallery modes both in circular 
and deformed microcavities are presented. In the case of a circular cavity, both analytical solutions 
and asymptotic formulas are derived. The coupling efficiencies at different coupling regimes for 
cylinder incoming wave are discussed, and the maximum efficiency is estimated for the practical 
Gaussian beam excitation. In the case of a deformed cavity, the coupling efficiency can be higher 
if the excitation beam can match the intrinsic emission well and the radiation loss can be tuned by 
adjusting the degree of deformation. Employing an abstract model of slightly deformed cavity, we 
found that the asymmetric and peak like line shapes instead of the Lorentz-shape dip are universal 
in transmission spectra due to multi-mode interference, and the coupling efficiency can not be 
estimated from the absolute depth of the dip. Our results provide guidelines for free space coupling 
in experiments, suggesting that the high-Q ARCs can be efficiently excited through free space which 
will stimulate further experiments and applications of WGMs based on free space coupling. 

PACS numbers: 42.79.-e, 42.79.Gn, 42.55.Sa 
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I. INTRODUCTION 



Whispering gallery modes (WGMs) were first ex- 
plained by Lord Rayleigh [l| in the case of sound propa- 
gation in St Paul's Cathedral circa in 1878. WGMs also 
exist in light waves which are almost perfectly guided 
round by optical total internal reflection in low loss res- 
onators 0. High quality (Q) factor excessing 10 10 has 
been achieved. Combining the ultra-small mode volume 
with the high-Q factor, the light-matter interaction can 
be greatly enhanced in WGM microcavities. Therefore, 
optical WGMs are gaining growing attentions in a wide 
range of fields including ultra-sensitive sensors Q, low 
threshold lasers Q , frequency combs § , cavity quantum 
electrodynamics (QED) [6], and quantum optomechanics 
@, [H . The WGMs also affect the scattering on spherical 
particles Q , which was found to play an important role in 
the glory phenomenon [ldj . In the last twenty years, the 
high-Q optical WGMs have been reported in various mi- 
crocavities, such as liquid droplet fllj . microspher e 1131 . 
microdisk [la ], microcylinder [lj] and microtoroid [la ]. 

The excitation and collection of the WGMs are es- 
sential issues in practical applications. It is taken for 
granted that free space coupling to WGMs is very inef- 
ficient since the radiation loss of circular shaped micro- 
cavity is isotropic and low, which brings the difficulty of 
energy transference between outside and WGMs. Near 
field couplers, such as prism fiber taper [T3, EU and 
perpendicular waveg uide[l9|, enable high efficient excita- 
tions and collections to microresonators. Therefore, these 
near field couplers are widely adapted in experiments. 
The underlying mechanism of coupling process between 
WGMs and couplers have been well understood [20T422T ] 
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and many phenomena in the spectrum are reported, such 
as the analogue of electromagnetically induced transpar- 
ent (EIT) [23l - l25j . asymmetric Fano line shape fl9l l2rj|- 
HH and ringing phenomena [2t| . However, there are some 
limitations in the near field couplers, such as large foot- 
print, stabilization requirement, and power limit (high 
excitation energy would cause strong nonlinear effects in 
fiber tapers). 

On another hand, the asymmetric resonant cavities 
(ARCs) have been demonstrated to be able to support 
high-Q WGMs and give highly directional emission 
|38| . The underlying principle and interesting chaotic ray 
dynamics in such an open billiard have been studied ex- 
tensively [39l - l43j . Especially, the unidirectional emission 
cavities have been successfully designed [4445 if . which 
enable high efficient collection of high-Q WGMs through 
free space. As the reversal of emission, the focused beam 
in free space can excite the WGMs, which have been stud- 
ied in experiments [52TI57I ) . However, the basic features 
and potentials of this free space beam coupling strategy 
are not fully explored. 

In this paper, we present a theoretical study on free 
space coupling to high-Q WGMs in both circular and 
deformed microcavities. First of all, the cylinder waves 
incident to a circular cavity is solved analytically and 
asymptotically in the cylinder coordinates, and the an- 
alytical equations are well consistent with the standard 
input-output formulation. Then, we consider the exci- 
tation of the circular cavity by Gaussian beam for prac- 
tical applications, the maximal efficiency of about 20% 
can be achieved under phase matching condition. Af- 
ter that, the excitation of deformed cavity by Gaussian 
beam is studied by an abstract model of multimode in- 
teraction. The result shows that the ARCs not only gives 
good matching parameters, but also can balance radia- 
tion and non-radiation losses, which can be well adapted 
in experiments for high efficient free space coupling. One 
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interesting result is that the asymmetric spectra is uni- 
versal in the free space coupling as a result of multimode 
interference, and the energy transferring can not be de- 
duced from the dip depth in the transmission spectrum. 



II. CIRCULAR CAVITY 

A. Cylinder wave 

We start with a two-dimensional circular shaped mi- 
crocavity, where electromagnetic field (tp) can be solved 
in the cylinder coordinates (r, <j)) analytically. Any 
electromagnetic field can be decomposed in the ba- 
sis of Bessel functions J m (kr) and Hankel functions 

(kr) with integral m 6 (—00,00). Since the 
field intensity is finite everywhere including the origin, 
thus the field inside the cavity should be represented 
by '^2 l a m J rn {nkr)e lm ^ when r < r c with r c is the radius 
of the cavity. The field outside (r > r c ) is represented 
by J2cmH^\kr)e im ^ and ^6 m ff r ( „ 1) (fcr)e"" < ^, which cor- 
respond to the inward and outward traveling cylinder 
waves, respectively. 

For m-th cylinder wave c m H$ {kr)e lm ^ incident to the 
cavity, only m-th components of cavity field and reflected 
outgoing wave are nonzero, due to the conservation of 
angular-momentum in this axial symmetry system. Ap- 
plying the boundary conditions, the coefficients of cavity 
field and outgoing field can be solved analytically with 



_ H£\z)H™\z) - h£>\z)H&\z) 
nPj; n (nz)H^\z) - J m (nz)H^>' (z) 
_ nPJ' m (nz)H£\z) - J m (nz)H^' \z) 
nPJ' m {nz)H^\z) - J m {nz)H^>'(z) 



(1) 

(2) 



where z = kr c , and p = —1(1) for transverse electric 
(TE) field (transverse magnetic (TM) field). 

In above formulas, there are singularities at z = zq 
where 

n p J' m {nz )H£> (z ) - J m (nz )H^'(z ) = 0. (3) 

This corresponds to the quasi-bound eigenmodes 
(WGMs) with dimensionless eigenfrequency (kr c = zo), 
which is a complex number. When the material refrac- 
tive index n is real, i.e. there no material absorption, we 
can write zq = Zq — ikq with Zq and kq are real numbers. 
Here, «o is the pure radiation loss of WGM, which cor- 
responds to the width of the resonance in spectrum as 
kq = Zq/2Qq with Qq is the radiation quality factor. 

For nearly resonant frequency z = zo + A (|A| <C \zq\), 
by expanding Eq. to the first order of A, we can 
approximately obtain 

nPJ' m {nz)H£Xkz) - J m [nz)H$\z) « F(z )A, (4) 



with F{z Q ) = (nP~ ] 
r<ipllj m (nz )H^(z )} 



- l)[(nJ^(nz )H^'(z ) + 
(nP +1 ^l)J m (nz Q )H^(z ). 
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FIG. 1: (color online) (a) The electric field distribution when 
a cylinder cavity is critical coupling with m — 30 cylinder 
wave, and detuning 8 — 0. (b) The normalized reflectivity 
and the total Q factor against the angular number m, with 
the refractive index n = 1.45 + 0.000117i. (c) and (d) are 
the intensity and phase of outgoing field for different cou- 
pling condition, with ki/kq = 0.2 (Red Solid lines), 1.0 (Blue 
Dashed lines) and 5.0 (Black Dotted lines). 



For the real materials, the non-radiation linear loss 
should be taken into account. For example, the material 
absorption can be considered in the characterize equa- 
tion by adding a small complex number (n^ <C n) to the 
refractive index as n = n + irii. For a monochromatic 
light with z ~ Zq — S, we have nz = (n + irii)(zQ — 5) ~ 
u(zq — S — ikq — ini) with k\ = zorii/n denotes the non- 
radiation loss. Substitute A = i(iS+KQ+Ki) into Eq. (0}, 
Eq. (TTJ), and Eq. ©. Under the high-Q mode (k < ^o) 
and near resonance (5 <C Zq) conditions, and with the ap- 
proximation that most of WGM energy is confined inside 
the cavity, we can get 

(X\ - 1 - 4 



b m (6) 



iS + kq + ki ttzqF(zq) 

iS - Kq + Ki F*(Zq) 



(5) 
(6) 



id + Kq + Ki F{Zq) 

On the other hand, the standard input-output formu- 
lation for a cavity mode is [58] 

d 



dt 



-E„ 



(-iS- 



Kq 



Ki)E n 



(7) 



where E m is the cavity mode field, E™ is the incoming 

(2) 

field in the form of radial inward cylinder waves Hm (zo), 
and the output field is £™* = -E™ + V2k^E m . When 



the system is in the steady state, 
satisfied, thus 
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should be 



(8) 



and the output field is 



Ei 



iS — Kq + Ki 
iS + Kq + Kl 



El 



(9) 



Comparing those equations with Eq. ([5]) and Eq. ([6]), 
the cavity and outgoing fields deduced by input-output 
formulation are consistent with the results derived by 
boundary conditions, only different in constants. The 
conversion relationships between two sets of formulas 



are dm — 



E': 



and 
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For the case of the near field coupler, the cavity field 
reads [H[ 




E m = 



id - 



(10) 



where Ki is intrinsic loss including radiation and absorp- 
tion losses, and n e is external loss induced by coupler. 
Comparing with Eq. ©, the free space and near field 
coupling manners are similar, with kq and K\ replaced 
by Ki and K e . 

From Eq. ((5|), the extremum of cavity field emerges 
when dE/dno = 0, i.e. Ko = n with 5 — 0. Similar to 
the case of a near field waveguide coupling to WGMs, 
there are three coupling regimes: under coupling regime 
Ko < K i; critical coupling regime kq = k\, and over cou- 
pling regime Ko > n\. The electric field distribution in 
Fig. 1(a) clearly shows the excitation of WGM at critical 

coupling by an inward cylinder wave (fcr)e l30( ^ with 
the spiral-like propagation. No reflection wave is excited 
in this case, which means all energy is perfectly absorbed 
by the cavity. Fig. 1(c) and (d) are the intensity and 
the phase of outgoing wave against the detuning in three 
regimes. In the under and over coupling regime, the en- 
ergy transfer efficiency without detuning is low. More- 
over, the phase of the outgoing wave is strongly changed 
by the WGM in the over coupling regime with small de- 
tuning value, while there is only small perturbation for 
under coupling regime. 

In the case of near field coupler, the extra coupling loss 
K e can be controlled by simply adjusting the gap between 
coupler and cavity, while Ki is an intrinsic parameter that 
does not change. Differently, in the free space coupling 
case, Ki is an intrinsic parameter determined by material 
and fabrication which is almost a constant for WGMs. It 
is only possible to change kq by change the cavity size 
or boundary shape. As shown in Fig. 1(b), the Q-factor 
[Q = kr c /2(K + Ki)] and the outgoing intensity change 
against m with n = 1.45 + 0.000117i. The Q-factor in- 
creases with the increasing m and shows a saturation 
value of about 6 x 10 3 . The critical coupling (kq ps k\) 
happens when to = 30, corresponding to a specific cavity 
size (r c ~ m/nk). For larger or smaller cavity size, the 
coupling efficiency reduces. 



B. Gaussian Beam 

Gaussian beam is most widely applied in experiments. 
Thus we use this beam in following studies for realistic 
consideration. Supposing a Gaussian beam centered at 



FIG. 2: (color online) (a) Schematic illustration of a Gaus- 
sian beam incident to a whispering gallery microcavity. The 
direction of beam is defined by 9q, and the location of center 
is (ro,(po) in the cylinder coordinator, (b) The electric field 
intensity distribution when the incident Gaussian beam is res- 
onant with the WGM, with kr c = 30.10, r = l.lr c , </f> = tt/2 
and 8q — 0. 



(xq, yo ) with intensity I — 1 and propagation along x-axis 
with the waist width w, the electric field at the waist is 



(v-vor 

G{x ,y) = e ^ 
Correspondingly, the angular spectrum is 



(11) 



— iyksin9 



e dy 



20? 



e 4 



(12) 



with e ( — f ) f )■ So, in the cylinder coordinator in Fig. 
2(a), the field distribution can be expressed as 



E(r, (/>) 



where 



W e -^ si ^e e ik( S in9y + co S 0x) d ^ ksine)j 

k 2V 7r 

(13) 



x = ?-cos(0 + 9 ) - r o cos(0 o + 9 ), (14) 
y = rsm((j) + 6 ) - r sin(</>o + #o)- (15) 

Substituting them to the Eq. (TT"3")) and employing the 
Jacobi- Anger expansion, 



zfcrcos((/>+(?o — 0) 



J m {kr)e 



im[-K/2-(<f,+9o-9)] 



, (16) 



m— — oo 



we can derive the coefficient of Gaussian beam in inward 
cylinder waves expansion 



im(ir/2-e )-ikq- [fc , ro ,f '" ] 



2^/1 + ^ 



(17) 



in the limit of kw 3> 1 and d = r o sin(0o + #o) an d 
q = rocos(</>o + ^o). Here, only the part of sin6* -C 1 
contributes in the integral, and the approximation is ap- 
plied 

cos(9 = Vl - sin 2 (9 w 1 - isin 2 6». (18) 

At 
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Substituting the coefficients to Eq. (1) and (2), the 
intracavity field and outgoing field can be solved ana- 
lytically for a cylinder cavity. For example, Fig. 2(b) 
shows the near field distribution when a Gaussian beam 
incident to the cavity with the on-resonance frequency to 
WGM with m = 40. Obviously, most energy are directly 
transmitted. Only a small portion of energy can be cou- 
pled into the cavity field. By integrating the Poynting 
vector of the electromagnetic field, we can calculate the 
power (energy flux) contained in the Gaussian beam 



Here, 



Pa = J- 



kw 



(19) 



2 2uj^ 

where u> is angular frequency, and /i is relative permeabil- 



1(2)) 



(kr)e 



im<fi 



ity. Power of normalized cylinder wave H., 
is 

Pm = — ■ (20) 

Therefore, the ratio of power contained in the m-th cylin- 
der waves to that in the Gaussian beam is 



\Cr, 



1 2 Pm 

1 iV 



(21) 



From Eq. (|17[) , rj m is maximal when q = and d = m/k. 
The first condition guarantees the smallest expansion of 
cylinder waves in momentum space (m) . The second con- 
dition corresponds to the phase matching condition that 
the momentum of the pump beam (kd) should be equal 
to the momentum of the WGM (m) . In practical applica- 
tions, we can manipulate the beam position d to satisfy 
the conditions, then 



Vn 



2 1 

7T kw 



(22) 



Therefore, the rate of power transfer from a Gaussian 
beam to the cavity WGM is limited to rj max . Under the 
paraxial approximation, the waist w of a Gaussian beam 
should be larger than 4//c, i.e. kw > 4, corresponding 
to the energy transferring rate r\ m < 0.2 for a cylinder 
cavity. 



III. DEFORMED CYLINDER 

A. General Model 

When the cavity boundary is slightly deformed (30l — 
[H, [45M50ll or small perturbations are introduced to the 
cavity [44l IHlj . highly directional emission can be ob- 
tained. Because the boundaries of ARCs are not regular, 
different momentum (m) components can couple to each 
other. In addition, WGMs always exist in pairs, clock- 
wise (A c ) and anti-clockwise (A a ) traveling modes, in a 
circular-like shape microcavity, then the general Hamil- 
tonian of ARC reads 



% — ft* 



Hi 



(23) 



M i( A L A 0,c + ^,a A j,a) (24) 



represents the energy of cavity modes, with Sj is the fre- 
quency difference between cavity mode and excitation 
laser, and Aj c (A,-. c ) is the creation (annihilation) opera- 
tor of cavity mode. The coupling between different modes 
is 

Hint = E y^X9i,k,cA\ c A ktC + gj,k,a,A] ia A k!a + h.C.) 

3 k>j 

+ J2Wj,kA{c A k,a + h.c), (25) 

j-k 

where f3j^ is the coupling coefficient between counter 
propagation WGMs and gj.k, a ( c ) IS * ne C0U P im g coef- 
ficient between co-propagation modes. The external 
pumping on modes is given by 

Hpump (^hj,rn,cAj c U m -\- ihj m a Aj a U_ m -\- /i.e.). 

(26) 

Here, hj^ mfi is the coupling strength of the cavity mode 
to the outside cylinder waves, with m — 1 ■ ■ ■ 00, where 
the external incoming field and outgoing field are rep- 
resented by cylinder waves as ^u™e* Vm Hm\kr)e~ lm ^ 



and J2 u ™ te ~ liPmH rr> V(kr)e- vm ' t > with e 
—Hm\kr c )/Hm\kr c ). Actually, A c and A a are 
time reversal to each other, so the coupling coefficients 
satisfy //,/... = g* k a = gj, k and % m , c = h* m a = h^ m . 

Here, we only concern about the ARCs with smooth 
boundary and small deformation which can support high- 
Q WGMs. Then we can omit the coupling between 
counter propagation WGMs (fii j = 0) and the direct 
scattering induced transition between incoming and out- 
going cylinder waves. Thus the Heisenberg equation of 
cavity field are written as 



—A. 

dt J ' c 



- iS ^ t 9 3 ,kA 



k>3 



k,c 



k<j 



A 



-XjA 



(27) 



where 



Xj 



= i5j 



%\o + Kj,i, Kjfi = J2 m \ h j,m\ /2 and 
Kj t i are intrinsic radiation loss and non-radiation loss, 
respectively. Supposing there is only one high-Q WGM 
with j = 1 near resonance of the excitation, and other 
modes are low-Q or largely detuned, we can adiabatic 
eliminate the fast varying terms by 4fAj <c = for j > 2, 
i.e. 



.4 



j,c\j>2 



Therefore, 



dt 



-i-^A hc 

Xj 



-X1A1 



OO h . 

m=l Xj 



00 

rn — 1 



(28) 



(29) 



5 



Here, Xl = Xi 
and S, 



n e + iS e with n e 

5a i i2 



«j,0+Kj,l 



j> 2 i^rpr i5i j i 
= — X)j>2 j^Jp IPij'I ) corresponding to the ex- 
tra loss and frequency shift induced by low-Q modes, 
and hi m = hi, m — iY^j>2 9i,j~%~~ are effective coupling 
coefficients to outside for high-Q WGM. Denoting the 
coupling efficiencies and the incoming field by vectors 
h = {hi^ m } and u — {u™} with m = 1 • • • oo, the cavity 
field of steady state becomes 



N /xi, 



(30) 



where £ = h ■ u/ h \u\ is the beam matching parameter. 
The output for arbitrary incoming field reads 

„out _ „,in , f 4 \ " \ " h j,m h j,m' m , 
u m - ~ u m + JmA\,c + 2^2^ ~ W ™' ' ' ^ 



j>2 m' 



Xj 



where /„ 



h* l m - iJ2j>2 Y~h* m corresponds to the 
radiation of A\ tC . 

From Cauchy-Schwartz inequality, |£| < 1, the maxi- 
mum of matching can be achieved only when 



h*/ 



(32) 



This corresponds to the optimal beam to excite the 

WGM E inc (r,<f>) oc E m C m ^i, m ^™ ) (^)e"™ , which 
is the reverse of the radiation of anti-clockwise WGM 
E r a ad {r,$) oc Y^ n hi, m e- i v^H^\kr)e im *, meaning that 
the optimal excitation of clockwise WGM can be achieved 
by the reverse of the radiation of anti-clockwise WGM. 
When on resonance, the cavity field reads A\ jC — 



\u\ /(ki.o + + K e)i and the energy transferring 



depends on both £ and 



1,0 



Ki,i + K e)- All the 



parameters mentioned above can not be solved analyti- 
cally because of the irregular boundary shape. Realistic 
details about the free space coupling should be solved by 
numerical simulation [62J. In the following, we analyti- 
cally solve the free space coupling to ARCs by an abstract 
model without the specific boundary shape. In this sim- 
plified model, we assume directional emission is in the 
form of Gaussian beam ignoring the specific boundary 
of a ARC. The general underlying mechanism and basic 
properties are revealed with several reasonable approxi- 
mations of a near circular boundary shape. 



B. Coupling Through Barrier Tunneling 

In a very slightly deformed ARC, the coupling between 
high-Q and low-Q modes are weak, and the directional 
emission is mainly due to the direct tunneling [59]. The 
center of the emission beam is departed from the bound- 
ary, as a result of barrier tunneling [60] . We assume that 
the parameters of emission Gaussian beam are d t = 1.2r c , 
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FIG. 3: (color online) (a) Field distribution of directional 
emission of anti-clockwise WGM in a near circular cavity (The 
field is shown in Logarithm scale, and the intracavity field is 
not shown since the exact boundary shape is not known in 
our abstract model) . The emission is in the form of Gaussian 
beam with dt = 1.2r c , wt ~ 0.1r c and direction is shown by 
arrows, (b) Field distribution when an on-resonance Gaus- 
sian beam coupling to the cavity in the over coupling regime 
with «i = 0.1«o, d = 1.2, w = 0.2r c and direction is shown 
by arrows, (c) Beam matching parameters against the d for 
different w. (d) The far field intensity of outgoing wave when 
a Gaussian beam incident to a cyiinder, for excitation fre- 
quency off-resonance and on-resonance at difference coupling 
conditions, (e), (f) and (g) are spectra detected at different 
angle (4>) in far field. 



Wt — 0.1r c and direction along x-axis, as shown in Fig. 
3(a). Thus, the coupling coefficients can be solved from 
the mode emission 



hi. 



f)oe 



(33) 



where f) = y^o/ 
ization J] |hi,m| : 



kw 
2 2 ' 



corresponding to the normal- 



2k . 



We discussed the coupling in detail using a Gaussian 
beam from left with 0q = and on-resonance to a high- 
Q mode of kr c = 31.1. Firstly, we calculated the beam 
matching parameter £ versus the beam position d for dif- 
ferent beam widths w, as shown in Fig. 3(c). The result 
shows that good match with £ > 0.8 can be easily real- 
ized with good fault tolerant to the imperfect focus or 
align of beam in experiments. In the near integrated sys- 
tem, we can neglect the coupling to other low-Q modes 
(k ' c + k\ c > K e ), then A c w ^/2ki, / (ki,o + 
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Similar to the cylinder case, the critical coupling with 
K i.o = K i,i guarantees the maximum energy transfer. 

With these approximations, we can calculate the far 
field distributions by substituting parameters into Eq. 
(30). The far field distributions at different coupling 
regime are shown in Fig. 3(d) with w = 0.2r c . At crit- 
ical coupling, the outgoing energy is great reduced but 
larger than as a result of imperfect matching. In under 
coupling, there is a small reduction of energy in all di- 
rections. However, in over coupling regime, the peak at 
4> = is split into two peaks, which is attributed to the 
interference between the transmitted beam and the emis- 
sion of the cavity mode, as a result of the strong phase 
shift at the over coupling regime. 

The spectra are also calculated (Figs. 3(e), 3(f) and 
3(g)) in different far field angles. It is revealed that the 
line shapes strongly depend on the positions of detectors. 
At <f> — 0°, the resonances shows regular Lorentz-shape 
dips similar to traditional near field coupler. In contrast 
with Fig. 1(c), similar dip depths are given in the over 
and critical coupling condition regimes, which indicates 
that the dip depth can not tell us the energy transfer 
rate in the free space coupling. Asymmetric line shapes 
(Fano-like and EIT-like) are observed at <f> = 15° and 
30°, as a result of multiple cylinder waves interference 
formalized in Eq. (|3"Tj) . 



C. Coupling Through Refraction 

When the degree of deformation is large, the direc- 
tional emission is mainly due to the radiation from low-Q 
modes. This process is also known as the dynamical tun- 
neling in the study of ARC, where rays in high-Q mode 
tunnel into chaotic sea and refracted out at some specific 
region in phase space [3J, |35[. The center of Gaussian 
beam should be adjusted on the cavity as a result of re- 
fraction [6l|. For the refraction dominated emission, the 
extra loss induced by low-Q mode is much larger than 
the intrinsic radiation loss ( i.e. K e 3> «o)> thus we can 
neglect the direct tunneling loss. Assuming the emis- 
sion is Gaussian beam with the parameters d r = 0.9r c , 
w r = 0.1r c and direction along x-axis (Fig. 4(a)), we 
have 
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FIG. 4: (color online) (a) Field distribution of directional 
emission of anti-clockwise WGM in a near circular cavity (The 
field is shown in normal scale, and the intracavity field is not 
shown since the exact boundary shape is not known in our 
abstract model). The emission is in the form of Gaussian 
beam with d r — 0.9r c , w r = 0.1r c and direction is shown 
by arrows, (b) The field distribution when an on-resonancc 
Gaussian beam coupling to the cavity in the over coupling 
regime with Ki = O.Iko, d — 0.9, w = 0.2r c and direction 
is shown by arrows, (c) The far field intensity of outgoing 
wave when a Gaussian beam is incident to the ARC with 
Gaussian beam emission, for off-resonance and on-resonance 
under different coupling conditions. 



Similar to the procedure in previous section, we solved 
the outgoing field when a Gaussian beam incident from 
left with 9q — 0, w — 0.2r c , d — 0.9r c and on-resonancc 
to a high-Q mode with kr c = 31.1. As shown in Figs. 
4(b) and 4(c), the maximal far field intensity is located 
away from 4> = as a result of the refraction of incident 
beam. The far field distribution is similar to the case 
of tunneling coupling, and the EIT or asymmetric line 
shapes are presented in the spectra (not shown here). 



Xj 



j>2 



f) e e 



/~ \kdf — ml 



(34) 



I Si 



with 



,\ z \h„ 
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In addition, the interactions 



7T k w 
2 2 



Here, 

2 Em>2 IXwiF 

between low-Q modes are not taken into consideration. 
We also assume the deformation just slightly influence 
the low-Q modes for simplicity, thus hj m = if j 7^ m. 
Therefore 



9l,r 



Xn 



f) e e 



-im7r/2- 



[kdj_ 



(35) 



IV. DISCUSSION 

According to the above analysis, the coupling to high- 
Q WGMs through free space can be efficient. For the 
cylinder wave, the largest coupling efficiency can be 
achieved under the phase matching condition, and the 
maximum is limited by waist width of a focused Gaus- 
sian beam, as the ratio of power contained in the cylinder 
wave is limited. For a deformed cavity, the efficiency can 
be higher since the pump Gaussian beam can be well 
matched to the directional emission beam. 

The energy transferred to the cavity is also determined 
by the coupling regime, i.e. the ratio of radiation loss to 



7 



non-radiation loss. The strongest coupling happens when 
k ps ki. The radiation loss kq decreases exponentially 
with the increasing of cavity size, but the non-radiation 
loss K\ is almost constant. Therefore, in a large cavity, 
kq is much greater than K\, which will give rise to very 
inefficient coupling. This limitation can be removed in 
the deformed cavity, as kq can be tuned by adjusting 
the deformation. Usually, a large degree of deformation 
can also lead to better directionality of cavity emission, 
thus better beam matching. Therefore, we can always 
tune the deformation degree to achieve the largest energy 
transfer rate, similar to the tuning of gap between near 
field coupler and microcavity. 

In contrast to near field coupler to microcavity, the 
asymmetric line sha pes are universal in free space cou- 
pling to WGMs [50, [621 . In more realistic cases, maybe 
multiple high-Q modes near to each other in spectrum 
are involved, and then the interactions between high- 
Q modes are enhanced by the low-Q modes mediated 
mode-mode interaction. From the view of ray dynamics 
in ARC, the high-Q modes usually correspond to regular 
period or quasi-period orbits. Due to dynamical tun- 
neling, the high-Q modes can couple to chaotic modes 
which correspond to chaotic ray trajectories. There- 
fore, the chaotic sea mediated the indirect coupling be- 
tween separated regular orbits in phase space (high-Q 
modes), which is also known as chaos assisted tunnel- 
ing [il|, [63|, H3] ■ Similar to the case of waveguide cou- 
pled microcavity [25l . [26j j , we can expect more profound 
asymmetric spectra when take more high-Q modes into 
consideration. 



V. CONCLUSION 

In summary, we present the theoretical study of the 
free space coupling to high-Q WGMs in both regular 



and deformed microcavities. The coupling by free space 
beam depends on the phase matching and beam match- 
ing conditions in regular cylinder cavities and ARCs, re- 
spectively. Three coupling regimes of free space coupling 
are discussed: the cavity should work near the critical 
coupling regime for high efficiency energy transfer. In 
cylinder cavity, the maximum energy transfer efficient is 
limited, and critical coupling can only be achieved with 
specific cavity size. The ARCs not only give good beam 
matching between focused Gaussian beam and highly 
directional emission, but also enable the tuning of the 
cavity degree of deformation to achieve critical coupling. 
Therefore, the efficiency approaching unity is possible in 
realistic experiments of free space coupling to unidirec- 
tional emission cavity. It is found that the asymmetric 
spectra or peak like spectra instead of the Lorentz-shape 
dip is universal in spectra, and the coupling efficiency 
cannot be estimated from the absolute depth of dip. Our 
results provide guidelines for free space coupling to high- 
Q WGMs, which will be valuable for further experiments 
and applications of WGMs based on free space coupling. 
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